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ABSTRACT 

This  report  outlines  an  investigation  of  an  ultrasonic  Acoustic  Electric  Feedthrough  (AEF)  to 
transfer  power  through  a  metal  plate.  This  approach  is  being  explored  as  a  potential  means  of 
wirelessly  powering  in  situ  structural  health  monitoring  systems  embedded  within  aircraft  and 
other  high  value  engineering  assets.  A  numerical  model  is  developed  describing  the  coupled 
electro-mechanical  behaviour  of  an  AEF  system,  and  is  then  validated  using  experimental  data. 
Usingapairofwell  matched  piezoelectric  disks  with  38  mm  diameter  and  2  mm  thickness,  and 
with  1W  of  real  input  power,  the  AEF  approach  is  shown  to  transfer  approximately  300  mW  of 
u  sabl  e  el  ectr  i  calpowerthroughaluminiumplatewiththickn  esses  intherangeofl.6mm-5mm. 
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Modelling  and  Experimental  Validation  of  the 
Acoustic  Electric  Feedthrough  Technique 


Executi  ve  Su  mmary 

ln-situ  Structural  Health  Monitoring  (SHM)  devices  allow  the  Australian  Defence  Force  to 
move  from  expensive  time-based  maintenance  approaches  for  ageing  platforms  to  cost- 
effectivecondition-based  approaches.  For  air  platforms  the  retrofitting  of  these  systems  is 
complicated  by  the  fact  that  the  majority  of  SH  M  devices  need  to  be  fitted  on  internal 
aircraft  structure,  underneath  the  aircraft's  skin.  IftheSHM  deviceisin  a  location  that  is 
difficult  to  access,  then  transferring  power  and  data  through  the  aircraft  skin  may  be 
necessary.  The  retrofitting  of  the  SH  M  device  then  becomes  complicated  due  to  flight 
worthiness  and  certification  issues,  with  thechief  concern  beingthat  hoi  esdrilled  through 
theaircraftskin  (to  passwiring)  could,  for  example,  act  as  stress  concentrators.  To  address 
this  issue  DSTO  has  investigated  the  use  of  an  ultrasonic  Acoustic  Electric  Feedthrough 
(AEF)  technique  to  pass  power  and  communications  through  metal  plates.  The  AEF 
techni que  consi sts  of  generati  ng,  transmitti  ng  and  recei vi  ng  hi  gh  frequency  (kH  z-M  H  z) 
acoustic  pulses  using  piezoelectric  elements.  A  numerical  AEF  modelling  capability  was 
developed,  and  used  to  explore  the  effect  of:  (i)  metal  plate  thickness,  (ii)  piezoelectric 
el  ement  thi  ckness,  (i  i  i )  the  type  of  pi  ezoel  ectri  c  materi  al ,  and  (i  v)  the  bond  I  i  ne  thi  ckness 
between  the  piezoelectric  element  and  metal  plate.  An  experimental  program  was  carried 
out  to  valid  ate  the  mod  el  predictions.  Experimental  studies  using  anon-optimized  system 
have  shown  that  the  AEF  approach  was  able  to  transfer  in  the  range  of  300  mW  of  real 
electrical  power  (for  1000  mW  of  real  input  power)  through  aluminium  plates  with  a 
thickness  range  of  1.6  mm  -  5  mm,  typical  for  aircraft  skin.  For  comparison,  a  standard 
recharge  on  an  80  mA-hr  lithium  polymer  battery  requires  four  hours  of  charging  at  a 
power  level  of  approximately  148  mW,  indicating  that  the  demonstrated  300  mW  AEF 
power  transfer  would  be  more  than  sufficient  to  recharge  a  smal  I  battery. 
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Nomenclature 


£ 33 ,rT 

undamped  relative  dielectric  constant  (dimensionless) 

£33S 

clamped  dielectric  constant  (F/  m) 

£0 

permittivity  of  free  space,  -8.8542x10 12  (F/  m) 

£ 33 

piezoelectric  constant  (C/  m2) 

kT,kp,  k33 

electromechanical  coupling  constants  (dimensionless) 

d31,  ds3 

piezoelectric  strain  constant  in  the  land  3  directions  respectively  {Cl  N) 

Y33E 

Young's  modulus  at  constant  field  (N/  m2) 

Cl3E 

closed  circuit  elastic  stiffness,  i.e.  electrodes  shorted  (N/  m2) 

C33E 

closed  circuit  elastic  stiffness,  (N/  m2) 

C33D 

open  circuit  elastic  stiffness  (N/  m2) 

si3E,  S33e,  sue  and  snr 

closed  circuitelasticcompliancesatconstantE  field  (m2/  N) 

Vd 

clamped  speed  of  sound  (m/  s) 
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mechanical  Q  (dimensionless) 
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density  (kg/  m3) 
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thickness  of  a  PZT  element  (m) 
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radius  of  a  PZT  element  (m) 

Ae 

area  of  a  PZT  element  (m2) 
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aspect  ratio,  height  to  radius,  of  a  piezoelectric  element 

f AR 

anti-resonant  frequency  (FH  z) 

fa 

resonant  frequency  (FI  z) 

AEF 

Acoustic  Electric  Feedthrough 

PZT 

Lead  ZirconateTitanate 
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1.  Background 


TheAustralian  Defence  Science  and  Technology  Organization  (DSTO)  is  developing  in  situ 
structural  health  monitoring  systems  [1,  2,  3,  4,  5,  6,  7,  8]  for  potential  use  in  high  value 
platforms  across  the  Australian  Defence  Force  (A  DF).  Using  smart  sensor  concepts,  airframe 
damage  and  damage  growth,  operational  strains,  corrosion,  corrosive  environments,  and 
other  structural  and  environmental  parameters,  would  becontinuously  monitored  on-board 
an  aircraft  to  provide  a  basis  for  real-time  damage  diagnostics.  This  technology  could 
potentially  permit  a  safe  reduction  in  inspection  and  regular  maintenancetimeand  costs  and 
thereforesubstantial  ly  i  impact  on  ai  rcraftthrough-l  ifesupport  costs  and  avai  labi  I  ity.  FI  owever 
providing  power  and  two  way  communications  to  a  smart  sensor  system  is  a  challenge 
particularly  when  physical  access  to  the  system  is  difficult.  To  address  this  issue,  DSTO  is 
i  nvesti  gati  ng  the  use  of  ul  trasou  nd  vi  a  an  A  cousti  c  Electri  c  F eedthrough  (A  E  F)  arrangement 
to  provide  power  and  communicati  onsthrough  theski  n  of  an  ai  rcraftto  smart  sensor  systems 
located  inside  the  aircraft. 

Hu  ct  nl  [9]  model  I  ed  an  A  EF  scheme  from  fi  rst  pri  nci  pi  es.  Sherri  1 6t  ul  [10]  at  the  NASA  Jet 
Propulsion  Laboratory  used  lumped  element  modelling  to  reproduce  Hu's  work,  and  to 
extend  the  modelling  of  the  AEF  arrangement  by  adding  a  more  realistic  description  of 
material  damping.  The  AEF  model  implemented  by  Hu  etal  [9]  and  Sherritefa/  [10]  consisted 
of  a  metal  plate  sand  w  i  ched  between  two  pi  ezoel  ectri  c  I  ayers.  Sherri  t  et «/  [10]  highli  ghted  the 
fact  that  their  modelling  did  not  incorporate  adhesive  bond  lines  between  the  piezoelectric 
material  andthemetal  pi  ate,  the  implicit  suggestion  being  that  the  bond  lines  might  affect  the 
behaviour  of  an  AEF  system.  This  report  examines  the  modelling  of  an  AEF  system  in 
LTSpice,  a  complete  and  fully  functional  SPICE  program  from  Linear  Technology  [11]. 
Implementing  the  model  in  LTSpice  allows  the  complexity  of  the  AEF  modelling  to  be 
su  bstanti  al  I  y  i  ncreased  by  al  I  ow  i  ng :  (i )  the  ad  d  i  ti  on  of  extra  A  E  F I  ayers  ( i  .e.  bond  I  i  nes)  and  ( i  i ) 
the  addition  of  realistic  el  ectri  cal  components  to  the  model.  Oncedevel  oped,  theLTSpiceAEF 
model  was  used  to  examine  various  geometrical  and  material  parameters  with  the  aim  of 
opti mizingthepower transferred  through  athin  metal  plate.Totestthefidelity  oftheLTSpice 
modelling  the  simulations  were  compared  to  experimental  data  collected  from  AEF 
arrangements  constructed  to  transfer  power  through  both  1.6mm  and  5mm  thick  aluminium 
plate. 

1.1  LTSpice  Modelling  of  a  Single  Piezoelectric  Disk 

Prior  to  modelling  an  AEF  scheme  it  was  necessary  to  ensure  that  LTSpice  was  capable  of 
model  I  i  ng  the  el  ectromechani  cal  behavi  our  of  a  si  ngl  e  pi  ezoel  ectri  c  d  i  sk.  There  are  nu  merous 
exampl  es  of  I  u  mped  el  ement  and  spi  ce  model  I  i  ng  of  pi  ezoel  ectri  c  el  ements  avai  I  abl  e  i  n  the 
literature[e.g.  10, 12, 13, 14, 15, 16].  Typically  these  are  based  on  either  the  Mason  [17]  or  the 
Krimholtz,  Leedom  and  M  atthaei  (KLM )  piezoelectric  models  [  18].  Both  theM  ason  and  KLM 
piezoelectric  models  are  one-dimensional  and  describe  the  thickness  mode  resonances  of  a 
piezoelectric  element,  utilizing  a  lumped  element  electrical  circuit  analog  to  model  the 
electromechanical  behaviour  of  the  el  ement. 
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The  LTSpi  ce  model  I  i  ng  of  a  Lead  Zi  rconateTitanate  (PZT)  el  ement,  as  descri  bed  i  n  thi  s  paper, 
requ  i  res  f  i  ve  p  i  ezoel  ectr  i  c  materi  al  parameters:  £33s,  e33,  c33D,  p  and  Q .  Some  cal  cu  I  ati  on  may  be 
required  to  implement  a  particular  material  typewithin  LTSpice  since  PZT  manufacturers 
rarely  specify  all  of  the  required  material  parameters.  The  physical  parameters  used  inthese 
model  s  vary  depend  i  ng  on  w  hether  the  pi  ezoel  ectri  c  el  ement  i  s  free  or  clamped.  A  n  exampl  e  of 
a  free  el  ement  is  one  hanging  freely  in  air.  Two  examples  given  in  the  literature  of  clamped 
elements  are:  (i)aPZT  el  ement  bonded  to  a  6  mm  thick  steel  plate  [10]  and  (ii)aPZT  element 
bonded  to  a  Poly  Methyl  Methacrylate  (PM  MA)  acoustic  load  matching  element  [15].  The 
parameters  are  also  dependent  on  geometry.  Kim  et  al  [19]  used  the  aspect  ratio  G  to 
investigate  the  effect  of  different  piezoelectric  disk  geometries:  e.g.  a  pancake  defined  by 
G  « 1 ,  and  cylinder  defi ned  by  G  » 1 .  Often,  key  parameters  I i kee33s,  e33,  and  c33D  will  need 
to  be  calculated  from  other  given  piezoelectric  material  parameters.  A  series  of  equations  that 
show  how  to  cal  cu  I  ate  these  parti  cu  I  ar  parameters  are  gi  ven  i  n  A  ppend  i  x  A .  These  equati  ons 
draw  heavily  on  the  various  piezoelectric  relationships  highlighted  by  Berlincourtef  al  [20]. 


The  LTSpice  model  of  a  PZT  element  shown  in  Figurelwas  adapted  from  the  literature  [10- 
16]  and  has  three  ports:  two  mechanical  ports  corresponding  to  the  back  (B)  and  front  (F)  of 
the  el  ement,  and  a  si  ngl  e  el  ectri  cal  port  (E)  with  the  assumpti  on  that  the  opposite  si  de  of  the 
PZT  element  is  grounded.  PZT  element  mechanical  losses  are  modelled  using  a  lossy 
transmission  line  as  outlined  in  the  literature  [14, 15].  The  LTSpice  library  and  schematic  files 
for  the  model  shown  in  Figure  lb  aregi  ven  in  Appendix  B,  with  parameter  definitions  given 
in  Appendix  C. 

1  direction 
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Figure  1:  (a)  Schematic  of  a  piezoelectric  disk  in  air  and  (b)  the  equivalent  LTSpice  model 
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Figure  la  shows  a  schematic  of  a  piezoelectric  disk  element  of  diameter  D  and  thickness  H 
w  here  the  resistors  RBT  and  RBR  model  the  acoustic  energy  lost  to  the  surrounding  air.  Rlis 
equal  to  50  Q  which  is  the  typical  output  impedance  of  an  AC  voltage  source.  Most  results 
presented  in  this  paper,  measured  and  modelled,  were  produced  using  piezoelectric  elements 
with  nominal  dimensions  of  D  =38  mm  and  H  =2  mm.  According  to  Kim  etal  [19],  for  PZT-5, 
if  the  aspect  rati o  G  ~ 0.1  there  shoul d  be  little  or  no  mode  coupl i  ng  between  the  radi al  and 
thi  ckness  vi  brati  on  modes  and  the  el  ectromechani  cal  cou pi  i  ng  constant  used  shou I  d  be  that 
for  a  pancake  geometry,  i.e.  kT  -0.47  [20].ThePZT  disks  used  in  the  current  work  had  an 
aspect  ratio  of  approximately  0.1,  yet  a  measured  [21]  impedance  spectrum  (Figure  2) 
revealed  substantial  modecoupling  between  thethi  ckness  and  radial  modes,  markedT  and  R 
respectively.  The  observed  mode  coupling  effect  will  bediscussed  later. 

ThePZT elementsused  inthispaperwereof  amaterial  composition  denoted  Pz27 [22],  which 
is  a  ceramic  with  material  properties  similar  to  PZT-5A.  The  model  ling  results  presented  in 
this  paper  were  calculated  using  a  Pz27  electromechanical  coupling  constant  of  k33  =0.699, 
where  the  33  subscript  denotes  a  cylindrical  geometry.  Interestingly,  mod  el  ling  the  element 
usi  ng  a  pancake  geometry  (w  i  th  kT  -0.47)  produced  a  su  bstanti  al  I  y  poorer  fit  to  the  measu  red 
data.  It  is  worth  noting  that  measured  values  of  modelling  parameters  were  used  where 
possi  bl  e.  A  I  i  st  of  parameters  used  for  the  pi  ezoel  ectri  c  cal  cul  ati  ons  i  s  gi  ven  inAppendixC. 


910 
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Figure  2:  Admittance  spectrum  measured  for  a  Pz27  disk  with  38  mm  diameter  and  2  mm  thickness 
showing  the  thickness  resonant  mode  T  and  the  coupled  radial  modes  R  [21] 
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1.2  LTSpice  Modelling  of  the  Acoustic  Electric  Feedthrough  Technique 

As  mentioned,  using  LTSpice  allowed  the  fidelity  of  theAEF  modelling  to  be  substantially 
i  ncreased  over  that  presented  i  n  the  I  iterature  by  al  I  owi  ng:  (i )  theadditi  on  of  extra  A  EF I  ayers 
(i.e.  bondlines)  and  (ii)  theadditi  on  of  realistic  electrical  components  in  themodel.  Figure3  is 
a  schematic  cross-section  of  an  AEF  system  showing  thetransmit  and  receive  PZT  elements 
bonded  either  si  de  of  a  metal  plate.  Figure4and  5show  schematics  of  two  AEF  modelsthat 
were  implemented  in  LTSpice.  Figure4  depicts  the  model  developed  to  extend  thework  of 
Flu  etal  [9]  and  Sherritef  al  [10].  Figure  5  shows  a  schematic  of  the  LTSpice  model  developed 
to  examine  acoustic  energy  transfer  through  an  aluminium  plate,  and  included  the  addition  of 
a  non-linear  electrical  load  on  the  receive  side  of  themodel.  Inclusion  of  realistic  bondline 
thi  cknesses  (i.e.  show  n  as  the  components  LossyBondTline  i  n  Fi  gure  4  and  5)  i  nto  the  model  I  i  ng 
was  required  to  achieve  reasonable  correlation  with  measured  AEF  power  transfer  results, 
which  will  be  discussed  later  in  this  report. 

To  perform  a  simulated  frequency  sweep  with  LTSpicea  script  was  written  to  execute  the 
following  simplealgorithm: 

•  modify  the  relevant  parameter  in  theLTSpicenetlist  (e.g.  frequency,  pi  ate  thickness, 
piezo  thickness), 

•  call  LTSpice  to  perform  a  transient  analysis, 

•  read  the  specified  LTSpice  measurements  from  the  LTSpice  log  file, 

•  record  the  measured  data, 

•  repeat  until  last  frequency. 

A  simulated  frequency  sweep  over  the  range  600  kFIz  -1300  kFIz  with  5  kFIz  steps  would 
typical  I y  take  24  hours  on  a  3  GFI  z  Penti  um  4  duo-core  machi  ne  with  3  Gb  of  fast  RA  M . 


receive 
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Figure  3:  Schematic  showing  an  AEF  configuration  in  cross-section 
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Figure  4:  LTSpice  model  of  a  Pz27/steel  AEF  system  with  similar  physical  configuration  to  that 
shown  in  Figure  3,  being  the  AEF  geometry  discussed  in  the  literature  [9,  10] 


Figure  5:  LTSpice  schematic  of  a  Pz27/Al  AEF  system  with  similar  physical  configuration  to  that 
shown  in  Figure  3  and  with  constant  1  W apparent  input  power  and  diode  bridge  RC  load 


Several  factors  were  found  to  havean  important  influence  on  thepower  transfer  efficiency  of 
the  AEF  system,  listed  below: 

(i)  Adhesive  bondline.  Asmentioned  earlier,  Sherritefa/  [10]  stated  that  the  adhesive  bond  line 
between  the  piezoelectric  element  and  the  metal  plate  might  limit  the  power  transfer 
efficiency  of  an  AEF  system.  In  this  report  AEF  systems  with  nominal  bondline  thickness  of 
100  pm  were  investigated  using  the  LTSpice  models  shown  in  Figure4and  5.  The  bond  lines 
areassumed  to  be  silver  loaded  epoxy  [23]  with  material  properties  given  by  Rajic  [24]. 

(ii)  Supply  impedance.  When  using  LTSpice  to  reproduce  the  steel  plate  case  reported  in  the 
literature  [9,  10]  it  was  assumed  that  the  supply  output  impedance  was  negligible  (wire 
resistance  only,  which  was  modelled  with  R2  =  0.1  Q  as  shown  in  Figure  4).  The  power 
ampl  ifier  utili  sed  experi  mental  ly  had  an  output  i  mpedance  of  50  Q  and  therefore  model  I  i  ng 
the  aluminium  pi  ate  AEF  system  required  R2  =  50  Q  (Figure  5). 

(iii)  Constant  power  source.  The  schematic  of  the  constant  power  sourcein  Figure5  is  shown  in 
Fi  gure  6  and  the  associ  ated  I  i  brary  fileis  included  inA  ppend  i  x  B.  The  constant  pow  er  source 
model  was  developed  to  allow  direct  comparison  between  experimental  observations  and 
model  predictions. 

Normally  1 W  of  apparent  input  power  was  chosen  to  drive  theAEF  systems  because  this 
could  easily  be  supplied  by  theamplifier  without  risk  of  damage.  When  thedrivefrequency 
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coi  nci  ded  w  ith  an  A  EF  system  resonance  the  phase  difference  between  the  i  nput  current  and 
voltage  was  near  zero  meaning  that  there  was  approximately  1  W  of  real  input  power  at 
resonance.  The  apparent  i  nput  power  was  measured  by  taki  ng  the  RM  S  of  the  i  nstantaneous 
product  of  the  i  nput  current  and  the  i  nput  voltage. 


{pDRJVE}  z1 
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>R4  -LT* 

1e8  50m 


Figure  6:  Schematic  of  the  constan  t  power  source 

(iv)  Material  damping.  Damping  in  the  piezoelectric  elements  is  modelled  using  a  lossy 
transmission  line  [14, 15].  The  LTSpice  models  used  in  this  report  extend  the  use  of  the  lossy 
transmission  line  to  the  metal  pi  ate  and  the  bond  lines. 

(v)  Diodes.  High  speed  Schottky  diodes  were  required  to  maximisetheAEF  power  transfer 
efficiency.  As  shown  in  Figure  5,  high  speed  IN  5818  diodes  [25]  were  used  to  create  the 
rectifying  bridge  on  the  receive  side  of  the  pi  ate.  The  IN  5818  Spice  model  supplied  by  the 
man  ufactu  rer  w  as  mod  i  f i  ed  to  i  ncl  u  d  e  measu  red  val  u  es  of  d  i  od  e  I  eakage. 

(vi)  Load  resistance.  To  duplicate  the  literature  results  the  investigation  of  a  steel  plate  AEF 
system  required  a  linear  20 Q.  load  resistor  as  shown  in  Figure  4.  A  30  Q  resistive  load  in 
parallel  with  a  high  frequency  1  pF  storage  capacitor  was  used  for  the  investigation  of 
al  umi  ni  um  pi  ate  A  EF  systems  (R1  and  C 1  respecti  vely  i  n  Fi  gure  5).  The  reason  for  choosi  ng  a 
30 Q  resistive  load  will  be  discussed  later. 

(vii)  Element  thickness,  H  =  2  mm.  Pz27  disks  with  38  mm  diameter  and  2  mm  thickness  were 
chosen  for  A  EF  i  nvesti  gati  ons  that  i  nvol  ved  al  umi  ni  um  plate.  The  2  mm  thi  ckness  was  chosen 
because  thefundamental  anti-resonant  frequency  was  approximately  1 M  Hz  which  was  low 
enough  to  ensure  that  cable-inductance  issues  were  negligible  and,  for  demonstration 
p  u  rposes,  si  mp  I  e  d  ri  ve  el  ectroni  cs  cou  I  d  befabri  cated  usi  ng  commerci  al  I  y  avai  I  abl  e  M  0  SFET 
driver  chips. 

1.3  Range  of  Validity  of  LTSpice  Modelling 

Thefol  I  owi  ng  secti  on  di  scusses  vari  ous  factors  and  restri  cti  ons  on  the  accuracy  of  model  I  i  ng 
an  AEF  systemin  LTSpice.  It  was  expected  thatthe  1-dimensional  LTSpice  models  presented 
would  produceuseful  predictions  as  long  as  the  total  acoustic  path  was  I  ess  than  theFresnel 
distance  for  the  solid  i.e.  the  AEF  system  operated  in  the  plane  wave  region.  Mason  [26] 
estimated  that  the  Fresnel  region  extends  for  a  distance, 
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For  theAEF  systems  investigated  in  this  paper,  with  a  drive  frequency  of  approximately  1 
M  FH z  and  u sing  a  Pz27  disk  with  a  diameter  of  38mm,  equation  1  predicts  a  Fresnel  distance 
of  d0  ~  90  mm.  Outside  of  the  Fresnel  region  diffraction  effects  need  to  be  included. 
Johansson  and  Martinsson  [27]  investigated  the  inclusion  of  diffraction  effects  in  PSpice 
simulations  of  ultrasonic  systems  using  an  approxi  mation  oftheequation  given  by  Kino  [28], 
This  describes  the  attenuation  between  two  identical  transducers  with  radius  R  and  separated 
by  d i  stance x  i  n  terms  of  the  Seki  [  29]  parameter  S, 


S 


(2) 


where  X  is  the  wavelength  of  the  ultrasound.  It  is  assumed  that  equation  2  holds  for  a 
continuous  wave.  Figure7plotstheattenuationasafunctionoftheSeki  parameter  and  shows 
that  for  S  <0.1  the  estimated  attenuation  due  to  diffraction  will  belessthan  0.59dB,  which 
represents  a  I  oss  of  13%.  For  the  cases  consi  dered  i  n  thi  s  report  the  Seki  parameter  i  s  I  ess  than 
0.1  indicating  small  diffraction  losses. 

The  radial  symmetry  of  a  piezoelectric  disk  allowed  Johansson  et  al  [27]  to  investigate  3- 
dimensional  diffraction  effects  with  a  1-dimensional  model,  however  it  is  more  difficult  to 
incorporate  other,  non-symmetric,  3-dimensional  eff  ects  i  nto  a  1-  dimensional  LTSpice  model. 
In  particular,  it  is  believed  that  small  non-uniform  variations  in  layer  thicknesses  (e.g. 
bondline,  piezoelectric  material,  metal  plate)  have  a  significant  effect  on  the  behaviour  of  an 
AEF  system  and  would  be  a  challenge  to  simulate  with  a  1-dimensional  model.  Due  to  time 
constrai  nts  average  thi  ckness  val  ues  have  been  used  i  n  the  current  model  I  i  ng. 


Figure  7:  Attenuation  caused  by  diffraction,  after  Johansson  [27]  and  Kino  [28] 

The  LTSp  i  ce  mod  el  I  i  ng  i  s  opti  mi  sed  f  or  the  f  i  rst  th  i  ckness  mod  e  as  d  escri  bed  by  P  uttmer  et  al 
[14],  Practi  cal  ly  thi  s  meant  that  a  constant  frequency-i  ndependent  resi  stanceval  ue  was  chosen 
for  each  of  the  lossy  transmission  lines  that  are  used  to  model  damping  in  the  piezoelectric 
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elements,  the  metal  plate,  and  the  bondlines.  Lossy  transmission  line  resistance,  RTl,  is 
calculated  using  [15], 


R 


TL 


®AR  L 

Q 


(3) 


where  L  =  pA,  A  =  xR2 ,  coAR  =2x  fAR  and  the  anti-resonant  frequency  fAR  =  vDl  (2  H)  as 
given  by  [20], 

The  assumption  of  a  constant  RTL  implied  by  equation  3  meant  the  mechanical  Q  values  used 
in  the  LTSpice  modelling  were  only  correct  at  the  anti-resonant  frequency  /AR,  and  the 
mechanical  Q  values  for  all  materials  effectively  increased  as  the  frequency  was  in  creased.  In 
other  words,  the  assumpti  on  of  a  constant  RTl  over-esti  mates  the  materi  al  dampi  ng  bel  ow/AR 
and  under-estimates  the  material  damping  above/AR. 

The  silver-epoxy  material  properties  determined  by  Raj  ic  [24]  were  calculated  at~4MHz.The 
LTSpice  model  ling  in  this  paper  assumed  that  the  mechanical  damping  of  silver  epoxy  was 
frequency  independent.  Theauthorsareawarehowever  that  thedamping  properties  of  visco- 
el  asti  c  materi  al  s  can  vary  quite  marked  I  y  w  i  th  frequency.  For  exampl  e,  van  Deventer  et  al  [  15] 
shows  that  the  mechani cal  Q  of  PM  M  A  i ncreases  from  10  to  100  over  the  10  H  z  - 10  M  H  z 
frequency  range. 

2.  Experimental 


Thi  s  secti  on  exami  nes  the  manufactureof  the  vari  ous  A  EF  arrangements  and  thesubsequent 
experimental  investigations.  Figure  8  shows  two  schematics  of  the  experimental 
configurations  utilised  to  examine  the  behaviour  of  the  A  EF  arrangements.  Two  well  matched 
Pz27diskswith  38  mm  diameter  and  2  mm  thickness  were  chosen  (wh  etrewell  matched  meant 
that  the  two  Pz27diskshad  well  matched  fundamental  thickness-mode  resonance  and  anti¬ 
resonance  peaks).  Silver-cement  was  used  as  the  electrode  material  for  the  Pz27  disks  because 
it  is  easily  soldered.  The  disks  were  bonded  to  either  side  of  the  metal  plates  using  silver 
loaded  epoxy  [23]  with  thedisk  centres  axially  collocated  to  within  +1  mm.  The  steel  plate 
used  had  dimensions  of  100  mm  by  64  mm  by  6  mm  thick,  whereas  the  aluminium  plates 
were  square  with  a  sidelength  of  600  mm. 

TheAEF  systems  were  prepared  using  thefollowing  procedure: 

(i) Impedance  measurements.  A  Solartron  Sl-1260impedanceanalyserwasused  with  acustom- 
built  PZT  measuring  attachment  (shown  in  Figure  9a  and  b).  Impedance  spectra  were 
measured  to  ensure  that  the  two  Pz27  disks  (Transmit  and  Receive)  were  well  matched  in 
terms  of  their  through-thickness  resonance  and  anti -resonance.  During  the  impedance 
measu  rements  the  sampl  es  were  mou  nted  w  i  th  the  I  oaded  gold  pogostick  contacti  ng  the  centre 
of  thedisk  (seeFigure9). 

Impedance  data  was  collected  over  the  frequency  range  lkFIz  to  lOMFIz,  sampling  at  1000 
points  per  decade,  with  1  second  temporal  averaging  per  point.  The  impedance  frequency 
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sweep  was  collected  near  the  thickness  resonance  at  maximum  frequency  resolution  to 
accurately  determinethefrequency  location  of  thethickness  resonance  and  anti -resonance.  If 
the  thickness  resonances  (and  anti-resonances)  of  the  two  disks  were  not  well  matched  (i.e. 
differed  by  more  than  +500  Hz)  then  the  process  was  repeated  until  two  matching  samples 
were  found.  Matching  of  the  two  disks  will  bediscussed  in  moredetail  later  in  this  report. 
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Figure  8:  Schematic  of  the  AEF  -power  transfer  experimental  arrangements  for  (a)  steel  plate  (b) 
aluminium  plate 


(ii)  Soldering  electrical  contacts  to  the  PZT  disks.  A  fine-point  sol  dering  iron  was  set  to  just  over 
360  °C .  FI  ux  was  appl  i  ed  to  a  smal  I  area  at  thecentre  of  a  Pz27  di  sk.  A  10  mm  I  ength  of  sol  der 
was  melted  onto  the  clean  soldering  iron  tip  and  thetip  was  then  quickly  'dabbed'  onto  the 
Pz27 surface.  Itwas  found  that  excessiveheat  or  del  ay  whilst  applying  thesolderwould  cause 
the  silver  metallization  (on  thePz27  surface)  to  be  subsumed  into  the  solder  ball. 

(iii)  Plate  surface  preparation.  The  area  of  plate  to  be  bonded  was  gently  cleaned  with  acetone 
and  followed  by  ethanol  to  removeany  acetoneresidue.  The  bonding  areawasthen  scored  in 
the  vertical  direction  using  600  grade  wet-and-dry  abrasive  paper  with  ethanol  lubricant 
(al  I  ow  i  ng  extra  area  for  the  si  I  ver-epoxy  spew  fi  I  let).  The  nu  mber  of  scori  ng  strokes  was  kept 
to  a  mi  ni  mum  to  reduce  materi  al  I  oss  from  the  surface  of  the  pi  ate.  Scori  ng  was  repeated  i  n 
the  horizontal  direction,  again  keeping  the  number  of  strokes  to  a  minimum.  Ethanol  was 
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used  to  remove  residuefrom  the  abrasive  paper.  The  bonding  area  was  cleaned  with  acetone 
and  again  the  acetone  residue  was  removed  with  ethanol. 

(iv)  Bonding  of  the  Pz27  disk.  The  CW2400  silver  loaded  epoxy  [23]  was  prepared  according  to 
the  manufacturers  instructions.  The  instructions  stated  that  the  pot-life  after  mixing  is  10 
minutes,  however  the  authors  found  that  the  useful  working  time  of  the  epoxy  was 
approximately  5  minutes.  Additionally,  the  instructions  state  that  mixing  should  take  2 
minutes,  however  for  small  quantities  the  authors  found  1  minute  to  be  sufficient,  which 
al  I  owed  extra  worki  ng  ti  me  for  pi  acement  of  the  Pz27  disk  onto  the  bondi  ng  area.  A  spatul  a 
w  as  u sed  to  ap  p I  y  a  u  n i  f or m  I  ay er  of  the  mi  xed  epoxy  to  th e  u  n d  ersi  d e  of  th e  Pz27  d  i  sk  w  h i  ch 
was  then  immediately  placed  on  the  plate.  Using  gloved  fingers,  pressure  was  applied  to  the 
disk  and  a  small  circular  motion  was  used  to  ensure  that  the  epoxy  spew  fillet  was  uniformly 
distributed  around  the  circumference  of  the  disk.  A  scalpel  was  used  to  remove  excess  epoxy 
f  rom  the  ed  ge  of  the  d  i  sk,  how  ever  the  spew-f  i  1 1  et  w  as  I  eft  i  ntact  d  u  ri  ng  th  i  s  p  rocess.  1 1  w  as 
found  that  removing  the  excess  at  this  stage  simplified  the  process  of  eliminating  electrical 
shorts  after  curing.  Figure  10  shows  an  example  of  a  38  mm  diameter,  2  mm  thick  Pz27disk 
bonded  to  an  aluminium  plate. 

(v)  Curing  of  the  silver  epoxy.  After  the  Pz27  disk  was  satisfactorily  positioned  it  was  covered 
with  Teflon  sheet  and  approximately  500  grams  of  mass  was  applied.  According  to  the 
manufacturer  the  epoxy  should  be  cured  at  room  temperature  (24  °C)  for  4  hours,  or  at  a 
temperature  of  66  °C  (to  121  °C)  for  10  minutes.  To  ensure  maximum  bond  strength  the 
authors  used  a250W  heat  lamp,  which  produced  a  sample  temperature  of  60°C,for  2  hours 
and  then  allowed  the  sample  to  sit  at  room  temperature  for  at  least  24  hours. 

(vi)  Removal  of  electrical  shorts.  After  curing  a  multi-meter  was  used  to  test  for  electrical 
shorting  between  thetop  Pz27  electrodeand  theplate.  If  a  short  had  occurred  a  Q-tip  and  a 
very  small  amount  of  acetone  was  used  to  remove  excess  epoxy  from  the  edge  of  the  Pz27 
d  i  sk,  w  i  th  care  taken  to  ensu  re  that  the  acetone  d  i  d  not  contact  or  f  I  ow  d  ow  n  to  the  spew  f  i  1 1  et 
and  hence  weaken  the  bond. 

Steps  (ii )  to  (vi)  were  repeated  to  bond  the  second  Pz27disk  to  the  opposite  side  of  the  plate. 
With  the  two  Pz27  disks  (transmit  and  receive)  bonded  to  the  plate  and  cured,  the  simple 
circuits  shown  schematically  in  Figure 4  (or  Figure  5)  were  built.  An  AG1021 100  W  power 
ampl  ifi  er  was  used  to  dri  vethetransmitti  ng  Pz27  di  sk.  A  4-channel  osci  1 1  oscope  with  fl  oati  ng 
probes  was  connected  to  monitor  the  input  voltage  and  current  used  to  drive  the  transmit 
Pz27  disk,  the  voltage  across  the  load  resistor  and,  if  a  diode  bridge  was  used  on  the  receive 
side,theoutputvoltageacrossthereceivePz27 disk  was  also  measured.  Drivefrequency  was 
step  p  ed  across  th  e  ran  ge  of  i  nterest  ( n  or  mal  I  y  400  k  FI  z  to  1300  khlz),withtheamplifieroutput 
adjusted  at  each  step  to  give  the  required  apparent  input  power  to  the  transmit  Pz27disk. 
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Figure  9: 


(a)  (b) 

Impedance  measurements  using  (a)  Solartron  1260  and  ( b )  the  DSTO  piezoelectric  disk 
adaptor 


Figure  10:  Pz27  element  immediately  after  silver-epoxy  bonding  to  aluminium  plate 
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3.  Results  and  Discussion 

In  the  following,  Part  1  develops  the  LTSpice  AEF  model,  extends  the  model  to  include 
bondlines  and  then  compares  the  predictions  with  those  found  in  the  literature  while  Part  2 
utilises  theLTSpiceAEF  model  to  investigate  geometrical  and  material  parameters  with  the 
aim  of  optimizing  power  transfer  through  an  aluminium  skin.  Part  2  also  includes  a 
comparison  with  experimental  data. 

3.1  Part  1  -  Steel  plate  Acoustic  Electric  Feedthrough 

TheAEF  geometry  under  investigation  is  that  considered  by  [9, 10]  and  involves  6  mm  thick 
steel  pi  ate  with  PZT  disks  of  1  mm  and  2  mm  thickness  bonded  to  either  side  of  the  pi  ate.  The 
LTSpice  model  in  Figure  4  was  developed  to  reproduce  this  and  to  also  examine  the  effect  of 
si  I  ver-epoxy  bond  I  i  nes.  M  ateri  al  parameters  for  mod  el  I  i  ng  w  ere  extracted  from  the  I  i  teratu  re 
and  are  given  in  Appendix  C.  The  LTSpice  model  was  optimised  for  the  thickness  anti¬ 
resonant  mode  of  the  1  mm  thick  transmit  piezoelectric  disk  located  at  a  frequency  of 
approximately  2  MFIz.  Figure  4a  shows  the  locations  of  V(out)  and  V(in)  and  Figure  11a 
compares  Vout  /  Vin  curves  reproduced  from  the  literature  [9]  with  those  produced  using  the 
model  given  in  Figure4. 


Figure  11:  Simulation  results  comparing  (a)  the  voltage  ratios  and  (b)  the  admittance  for  three  cases: 

(i)  literature  reproduced  from  Hu  et  al  [9],  (ii)  LTSpice  modelling  using  lossy  transmission 
line  damping  models  for  the  piezoelectric  and  steel  elements  and  (Hi)  the  LTSpice  model 
with  100  jum  silver  epoxy  bondlines  added  between  the  piezoelectric  elements  and  the  steel 


The  lossy  transmission  lines  utilised  in  theLTSpiceAEF  modelling  were  more  sophisticated 
than  the  damping  models  utilised  in  the  literature  [9,  10].  The  use  of  more  sophisticated 
damping  models  resulted  in  smaller  resonant  peaks,  as  shown  in  Figure  11a.  Figure  11a  also 
showsthe  effect  of  adding  two  100  pm  bondlines  (i.e.  one  bondline  on  each  side  of  the  steel 
plate)  to  the  LTSpice  model,  resulting  infurther  substantial  suppression  of  the  resonant  peaks 
(with  peak  voltage  ratios  VOUT  / VIN  <  I )  which  were  also  shifted  in  frequency.  Figure  lib 
showsa  comparison  of  theinput  admittance  curve  from  the  literature  and  those  modelled 
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usi  ng  LTSpi  ce  and  agai  n  the  additi  on  of  the  two  100  pm  bondl  i  nes  has  a  si  gnifi  cant  effect  on 
the  behavi  our  of  the  A  EF  model . 

Theinitial  AEF  modelling  was  carried  out  using  piezoelectric  material  parameters  extracted 
from  the  literature  [9,  10].  To  extend  the  modelling  and  hence  allow  comparison  with 
measured  AEF  data  it  was  necessary  to  show  that  the  LTSpice  piezoelectric  model  could 
accurately  predict  the  thickness  mode  resonance  behaviour  of  the  Pz27  disks  used.  The 
piezoelectric  LTSpice  model  shown  in  Figure  lbwas  hence  created  to  confirm  that  the  three 
port  p  isoelectric  mod  el  found  in  the  literature  was  suitable.  Notethat  the  LTSpice  model  was 
optimised  for  the  first  thickness  mode,  as  discussed  in  [14]. 

Impedance  spectra  from  the  model  were  compared  with  measured  impedance  spectra  that 
were  collected  using  a  custom  built  sample  holder,  shown  in  Figured  Figure  12  shows  the 
measured  impedance  data  from  a  free  Pz27disk,  38  mm  diameter  and  2mm  thick,  over  the 
range  1  kFH  z  to  10  M  H  z;  pi  otted  near  1 M  H  z  i  s  the  model  I  ed  i  mpedancedata  computed  usi  ng 
the  LTSpice  model  of  aPz27disk  shown  in  Figure  lb. 


frequency  (kHz) 


Figure  12:  Measured  impedance  magnitude  spectrum  over  the  frequency  range  1  kHz-10  MHz,  and 
modelled  data  in  the  vicinity  of  the  fundamental  through-thickness  resonance  (-972  kHz) 

The  series  of  measured  resonance  peaks  located  between  ~50  kHz  and  600  kHz  are  lateral 
vibrations  (in  the  1  direction)  that  are  not  modelled  in  the  LTSpice  Pz27  simulation,  which 
only  models  thickness- mode  resonances.  Figure  13  examines  the  frequency  range  near  the 
thickness- mode  resonance  and  shows  that  excellent  agreement  was  found  between  the 
measured  and  modelled  magnitude  and  phase.  Figure  13  also  shows  that  mode  coupling 
occurs  between  the  radial  and  thickness  modes.  N  oting  that  impedance  is  the  reciprocal  of 
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admittance,  to  fad  I  itatecompari son  of  Figure2with  Figure  13a,  thethicknessmoderesonance 
has  been  marked  withT  andtheradial  coupled  modes  are  label  led  R.  Asdiscussed  earlier,  the 
observed  mode  coupling  was  confirmed  by  the  manufacturer  [21]. 

H  avi  ng  confi  rmed  the  LTSpi  ce  pi  ezoel  ectri  c  model  was  pred  i  cti  vethe  model  was  then  appl  i  ed 
totheAEF  geometry  considered  in  [9]  and  [10],  after  which  an  experimental  validation  was 
performed.  Two  Pz27  disks,  38  mm  diameter,  with  thicknesses  of  1  mm  (transmit)  and  2  mm 
(receive)  were  bonded  to  either  si  de  of  a  smal  I  steel  pi  ate,  with  di  mensi  ons  100  mm  x  64  mm  x 
6  mm  thick.  A  20  Q  linear  load  resistance  was  attached  to  the  receive  Pz27  disk  as  shown 
schematically  in  Figure  4.  The  apparent  input  power  for  the  experiment  was  limited  to 
500  mW  as  higher  input  powers  were  found  to  distort  the  output  load  voltage. 


(a) 


(b) 


Figure  13:  Comparison  of  measured  and  modelled  (a)  impedance  magnitude  and  (b)  phase  near  the 
thickness  resonance,  where  thickness  mode  resonance  has  been  labelled  with  T  and  the  radial 
coupled  modes  are  labelled  R 
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Figure  14:  Oscilloscope  screen  capture  showing  three  traces  from  steel  plate  AEF  experiment.  Channel 
locations  are  indicated  in  Figure  8a:  Chl=V(in),  Ch2=l(R2),  Ch3=V(out)  across  20  Q 

Figure  14  shows  an  example  oscilloscope  trace  measured  at  a  drive  frequency  of  1075  khlz 
near  the  anti-resonance  frequency  of  the2  mm  thick  Pz27  disk  in  the  steel  plateAEF.  An 
oscilloscope  trace  was  captured  at  each  drive  frequency  step  shown  in  Figure  15.  Using  the 
oscilloscope  traces  itwas  possible  to  find  the  power  transfer  efficiency  v,  defined  as, 
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V  =  lOOx 


1  RECEIVE 

\Ptransmit  J 


(4) 


where  PTRANSMIT  V(///)x  f(R2)  end  P^ceive  V (out)  x  I  (R\) 

Figure  15  shows  a  comparison  of  measured  power  transfer  efficiencieswiththosefound  using 
LTSpice  simulation.  For  the  purpose  of  modelling,  nominal  100  pm  silver-epoxy  bond  lines 
were  chosen.  The  experimental  efficiency  was  calculated  using  PTRANSMIT  =500 mW 

and  PRECEIVE  =  )2  / 200 .  Figure  15  shows  that  both  the  measured  and  modelled  peak 

power  transfer  efficiencies  were  quite  poor,  less  than  2%,  for  this  geometry  and  load 
resistance. 

Subsequent  micrometer  measurements  revealed  that  the  thickness  of  the  steel  pi  ate  varied 
from  5.80  mm  at  the  edge  of  the  pi  ate  to  5.99  mm  at  the  mi  ddl  e.  Given  that  there  was  a  200  pm 
thickness  variation  across  the  pi  ate  it  was  not  expected  thatthefrequencies  of  the  measured 
resonance  peaks  would  match  those  computed  by  the  one-dimensional  LTSpice  simulation, 
which  assumed  uniform  bondlines  and  plate  thicknesses.  The  authors  expected  that  the 
thickness  variation  ofthesteel  platewould  mean  thatafull  three-dimensional  model  would 
be  required  to  satisfactorily  model  the  AEF  system,  however  Figure  15  does  seem  to  show 
reasonable  correlation  between  the  measured  and  modelled  resonant  peaks  at  625  khlz  and 
1025  kHz. 


2  | - 1 - 1 - 1 - 1 - 1 - r 


*  measured 
- LTSpice  simulated 


500  1000  1500  2000  2500 

frequency  (kHz) 


Figure  15:  Measured  versus  modelled  power  transfer  efficiency  for  a  1  mm  Pz27/6  mm  steel/2  mm 
Pz27  AEF  system.  Note  the  peaky  nature  of  the  plot  with  very  low  power  transfer  between 
the  peaks. 
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The  AEF  efficiency  curves  given  in  the  literature  [9, 10]  show  significant  power  transfer  in 
between  the  resonant  peaks,  however  with  more  realistic  damping  in  the  AEF  model  the 
efficiency  curves  became  peaky  with  very  little  power  transfer  between  theresonant  peaks.  It 
isal  so  worth  noting  that  the  efficiencies  shown  in  Figure  15  are  substantially  lower  than  the 
predictions  found  in  the  literature  so  it  is  clear  that  using  lossy  transmission  lines  to  model 
damping  and  adding  the  adhesive  bondlines  has  a  significant  effect  on  the  power  transfer 
efficiency. 

3.2  Part  2  -  Aluminium  plate  Acoustic  Electric  Feedthrough 

The  u I  ti  mate  goal  of  th i  s  research  i  s  the  d evel  op ment  of  an  A  E F  system  to  p rovi  d e  pow er  and 
communications  through  the  skin  of  an  aircraft.  With  the  LTSpice  AEF  model  agreeing 
satisfactorily  with  theliteratureand  alsoshowntohavepredictivecapability,  an  investigation 
of  AEF  power  transfer  through  aluminium  plate  was  undertaken. 

Prior  to  modelling  the  aluminium  plate  AEF  system  a  series  of  impedance  curves  were 
measured  beforeand  after  bonding  Pz27disks  (38 mm  diameter  and  2mmthick)toal.6mm 
thick  aluminium  plate.  Theresults  of  these  measurements  areshown  in  Figure  16.  Figure  16a 
shows  impedance  curves  taken  overthefrequency  range  lkhlz  to  10  M  FI  z  and  compares  the 
impedance  of  the  transmit  Pz27disk  beforeand  after  bonding  to  a  1.6  mm  thick  aluminium 
plate.  As  might  be  expected,  bonding  the  Pz27  disks  to  a  metal  plate  resulted  in  the 
attenuation  and  shifting  of  the  radial  resonance  peaks  located  between  50kFlzand  600kFlz. 
Thethi  ckness  resonance  modes  seen  at  1 M  FI  z  and  above  were  al  so  attenuated  and  shifted  (to 
a  higher  frequency)  after  bonding.  Figure  16b  showsthatthefundamental  thickness  resonance 
of  both  the  transmit  and  receive  Pz27  disks  shifted  to  a  higher  frequency  after  bonding,  i.e. 
from  972  kFH  z  to  -1060  kFH  z. 

C I  ose  exami  nati  on  of  Fi  gure  16b  reveal  s  that  the  anti -resonant  frequenci  es  of  both  thetransmit 
and  recei  ve  Pz27  disks  was  1084  kFIz  prior  to  bonding  and  shifted  to  1129  kFIz  after  bonding. 
Figure  16c  shows  that,  prior  to  bonding,  thetransmit  and  receive  Pz27  disks  had  extremely 
well  matched  resonance  and  anti -resonance  frequencies  with  the  centre  of  the  phase  peak 
located  at  -1027  kFH  z  and  that,  after  bondi  ng,  the  peak  shifted  to  -1100  kFH  z.  Pri  or  to  bondi  ng 
the  mode  coupling  between  thethi  ckness  and  radial  modesisdearly  visiblein  Figure  16b  and 
Figure  16c  which  also  show  that  the  mode  coupling  was  severely  attenuated  after  bonding. 

Simulations  werecarried  out  using  the  LTSpice  AEF  model  in  Figure5and  comparison  was 
madewith  laboratory  data  measured  using  theexperi mental  setup  (Figure8b).  Anexampleof 
a  measu red  osci  1 1  oscope  trace  i  s  show  n  i  n  Fi  gu  re  17;  si  mi  I  ar  traces  w ere  captu  red  at  each  step 
in  drivefrequency.  The  measured  data  presented  in  Figure  18, 19  and  20  was  calculated  from 
traces  similar  to  Figure  17. 
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(a) 


resonance 

(b)  (c) 

Figure  16:  Impedance  measurements  of  2  mm  thick  Pz27  disks  before  and  after  bonding  to  1.6  mm  Al 
plate:  (a)  shows  the  measured  impedance  magnitude  from  the  transmit  piezoelectric 
element,  for  both  transmit  and  receive  piezoelectric  elements  (b)  shows  details  of  the 
measured  impedance  magnitude  near  the  fundamental  thickness  resonance  and  (c)  the 
corresponding  phase 
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Figure  1 7:  Oscilloscope  screen  capture  showing  four  example  voltage  traces  from  the  aluminium  plate 
AEF  experiment.  Channel  locations  are  indicated  in  Figure  8b:  Chl=V(in),  Ch2=I(R2), 
Ch3=V(out)  across  30  Q,  Ch4=voltage  across  the  receive  piezo 

Figure  18a  is  the  model  led  and  measured  load  voltage  as  a  function  of  frequency  for  an  AEF 
formed  using  1.6mm  thick  aluminium  plate,  showing  a  modelled  peak  voltage  of  ~4  V  and  a 
measured  peakvoltageof  ~3V.  Plotted  in  Figure  18b  arecurves  of  load  power  asafuncti on  of 
frequency.  Si  nee  the  I  oad  voltage  had  a  negl  i  gi  bl  e  amount  of  ri  ppl  ethe  power  was  esti  mated 

usingthesimplerelationship  PL0AD  =  {v^an  )2  /30Q  .Themodelled  peak  power  transfer  was 
near  550  mW  whilst  the  measured  peak  power  transfer  was  just  over  300  mW.  Itisassumed 
thatthedifferenceisthepowerlosttothemetal  platewhich  isnotincluded  in  the  mod  el.  The 
approximate  power  transfer  efficiency  v  =  1 00  x  (PLOAD  / 1 ) ,  show  n  as  a  percentage  i  n  Figure 
18b,  is  only  accurate  at  the  resonant  peaks  w  here  the  real  input  power  was  about  1 W. 

Figure  18b  also  shows  that  the  maximum  measured  power  transfer  peak  was  centred  near 
1100  kFIz  with  a  Full-Width-Flalf-Maximum  (FWFIM)  of  -166  kFIz  while  the  maximum 
modelled  power  transfer  peak  was  centred  at  1091  kFIz  with  FWFIM  -183  kFIz.  Both  the 
modelled  and  measured  power  curves  show  that  the  centre  of  the  maxi  mum  power  transfer 
peak  was  located  near  the  anti -resonant  frequency  of  a  bonded  2  mm  thick  Pz27disk  (Figure 
16)  which  isexpected  si  nee  losses  within  the  disk  areminimised  at  anti-resonance.  The  phase 
pi  ot  between  input  current  and  drive  voltage  shown  in  Figure  19  gives  an  indication  astowhy 
the  resonance  peak  near  1100  kFIz  was  so  wide.  Themodelled  phasedata  shows  threecl  early 
discernable  peaks  located  side  by  side  near  1100  kFIz  which  indicates  that  there  may  be  a 
similar  number  of  peaks  present  between  lOOOkFIzand  1200  kFIz  in  the  measured  phasedata. 

Modelled  and  measured  results  pertaining  to  a  5  mm  thick  aluminium  plate  are  shown  in 
Figure  20  and  are  quite  similar  to  those  in  Figure  18  for  the  1.6  mm  thick  plate.  In  particular, 
the  peak  values  of  the  modelled  and  measured  load  voltage  (and  power)  at  -1040  kFIz  are 
close  to  those  found  for  the  1.6  mm  plate  near  1100  kFIz. 
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frequency  (kHz) 


(a) 


Figure  18:  Measured  and  modelled  results  for  the  2  mm  Pz27/1.6  mm  Al/2  mm  Pz27  AEF 
configuration  (a)  load  voltage  and  (b)  transmitted  power  as  a  function  of  frequency 


To  examine  the  effect  of  metal  plate  thickness  on  the  power  transfer  efficiency  a  series  of 
simulations,  varying  the  inputdrivefrequency  and  thealuminium  plate  thickness,  wererun 
using  the  LTSpice  model  shown  in  Figure  5.  Again,  38  mm  diameter,  2mm  thick  Pz27disks 
were  used  to  form  the  AEF  system.  Figure  21  shows  the  load  voltage  response  for  a  range  of 
fixed  drive  frequencies  for  varying  aluminium  plate  thickness.  From  Figure  21  it  can  be 
inferred  that  for  drive  frequencies  near  the  1100  kHz  anti -resonance,  the  thickness  of  the 
aluminium  plate  has  little  effect  on  the  amount  of  power  transferred  through  the  plate.  This 
prediction  isborneoutin  practice  as  shown  in  Figure  18  and  Figure  20  where  the  measured 
power  transfer  effi  ci  ency  was  found  to  be  si  mi  I  ar  for  al  umi  ni  um  pi  ate  thi  cknesses  of  1.6  mm 
and  5  mm.  Of  course  as  plate  thickness  increases  diffraction  effects  (i.e.  Figure  7)  become 
important  and  would  reduce  the  power  transfer  efficiency. 


frequency  (kHz) 


Figure  19:  Phase  spectra  for  a  2  mm  Pz27/1.6  mm  Al/2  mm  Pz27  AEF  configuration 
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Figure  20:  Measured  and  modelled  results  for  the  2  mm  Pz27/5  mm  Al/2  mm  Pz27  AEF 
configuration,  (a)  load  voltage  and  (b)  transmitted  power  as  a  function  of  frequency 

For  dri  vefrequenci  es  away  from  1100  kH  z  (i  .e.  for  frequencies  away  from  the  anti -resonance) 
thereappeared  to  bea  particular  set  of  thicknesses  that  allowed  maximum  power  transfer.  For 
example,  Figure  21  indicates  that  a  drive  frequency  of  900  kHz  produces  maximum  power 
transfer  through  plates  of  0.6  mm  and  3.5  mm  thickness. 


—i — i — i — i — | — i — i — i— 


—i — i — i — i — | — i — i — i— 


ftvy/ik  m 

V.  •  aa.4*  a. 


2  3  4 

Aluminium  thickness  (mm) 


Figure  21:  Modelled  results  for  an  AEF  system  with  varying  plate  thickness  and  drive  frequency  using 
alW  constant  apparent  input  power  and  a  30  Dload  resistor 

The  LTSpice  model  in  Figure  5  was  modified  to  investigate  the  effect  of  the  piezoelectric 
element  thickness  on  the  power  transferred  through  a  1.6  mm  aluminium  plate.  The  load 
resistor  R1  was  removed  and  the  voltage  across  the  1  pF  storage  capacitor  Cl  was  measured 
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as  a  function  of  time.  Both  the  transmit  and  receive  piezoelectric  elements  used  intheAEF 
system  had  the  same  nominal  thickness.  The  input  drivefrequency  was  adjusted  to  be  equal 
to  the  resonance  frequency  for  each  element  thickness  modelled.  Figure  22a  shows  the 
modelled  steady  state  load  voltage  as  a  function  of  Pz27  disk  thickness.  The  actual  voltages 
plotted  in  Figure  22a  are  90%  of  the  final  steady  state  Cl  capacitor  voltage,  V90%)  this  value 
was  chosen  to  keep  simulation  times  reasonable. 


piezo  element  thickness  (mm) 

(a) 


piezo  element  thickness  (mm) 

(b) 


piezo  element  thickness  (mm) 
(0 


Figure  22:  Modelled  results  for  an  AEF  system  for  a  1.6  mm  thick  Al  -plate,  where  the  PZT  element 
thickness  is  varied:  (a)  90%  of  steady  state  load  voltage,  (b)  time  taken  to  reach  90%  of 
steady  state  load  voltage,  and  (c)  load  power  as  a  function  of  varying  Pz27  thickness 


The  time  taken  for  the  Cl  capacitor  to  reach  90%  of  its  steady  state  voltage,  t90%,  is  plotted  in 
Fi  gure  22b.  Fi  gure  22c  pi  ots  I  oad  power  as  a  functi  on  of  pi  ezo  el  ement  thi  ckness,  where  I  oad 
power  was  calculated  using, 
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n  V  90%  ) 

Pload  =  -  (5) 

*90% 

whereC=Cl=l  pF. 

Theload  power  plotted  in  Figure  22c  is  not  power  in  thetruesense,  but  is  the  power  stored 
and  consumed  in  the  foil  owing  cycle: 

(i)  a  capacitor  charge  cycle  (where  the  load  is  capacitive  only)  is  followed  by, 

(i  i )  an  energy  use  cycle  duri  ng  whi  ch  the  energy  stored  i  n  the  capacitor  i  n  stage  (i )  i  s  di  ssi  pated 
through  a  real  load. 

Figure  22c  shows  a  trend  towards  lower  power-transfer  efficiency  for  piezoelectric  element 
thicknesses  below  1  mm.  For  a  1.6  mm  aluminium  plate  Figure  22c  suggests  that  a  range  of 
piezoelectric  element  thicknesses  from  1mm  to  3mm  could  beused  to  form  an  AEF  system 
with  power  transfer  efficiencies  in  the  range  of  25-40%. 


(a)  (b) 

Figure  23:  Modelled  results  for  2  mm  Pz27/1.6  mm  Al/2  mm  Pz27AEF  configuration,  (a)  load  voltage 
and  (b)  load  power,  with  varying  load  resistance  and  drive  frequency  swept  across  the 
fundamental  thickness  resonance,  1025-1150  kHz 


Finally,  the  reason  for  choosing  a  30  Q.  load  resistance  for  the  aluminium  plate  AEF 
configuration  is  apparent  in  Figure 23  which  plots  a  series  of  simulation  results  where  both 
theload  resistance  and  drive  frequency  were  varied.  The  LTSpice  model  shown  in  Figure  5 
was  again  modified  and  this  time  resistor  R1  was  stepped  across  the  range  IQ  to  100  Q  and 
the  input  drive  frequencies  were  stepped  from  1025  kFIzto  1150  khlz,  corresponding  to  the 
main  power-transfer  peak  as  shown  in  Figure  18  and  Figure  20.  Figure  23  shows  that  a  30D. 
load  resistor  produced  load  voltages  from  3.7  V  to  4.1  V  which  was  deemed  to  be  a  useful 
voltage  range  by  the  authors.  Flowever,  more  importantly  it  appears  that  the  load  resistance 
for  optimum  power  transfer  efficiency  lies  somewhere  in  the  range  of  30  Q  to  60  Q  and  is 
somewhat  i  ndependent  of  thedri  vefrequency  as  long  as  it  is  withi  n  approxi  mately  +50  kH  z  of 
the  thi  ckness-mode  resonance. 
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4.  Conclusion 

The  Acoustic  Electric  Feedthrough  (AEF)  technique  was  investigated  as  a  possible  means  of 
passi  ng  power  through  theal  umi  ni  um  ski  n  of  an  ai  rcraft.  A  n  A  EF  operates  usi  ng  two  axi  al  ly 
aligned  piezoelectric  elements,  located  on  opposite  sides  of  a  metal  plate.  A  piezoelectric 
el  ement  exci  ted  at  i  ts  th  i  ckn  ess  mod  e  anti  -  reson  ant  f  requ  en  cy  i  s  u  sed  to  tr an  smi  t  u  I  trasou  n  d 
through  a  metallic  plate.  The  ultrasound  is  received  by  a  second  element  located  on  the 
opposite  side  of  the  pi  ate.  An  LTSpice  model  oftheAEF  scheme  was  developed  and  used  to 
reproduce  and  extend  modelling  found  in  the  literature.  In  particular,  the  model  included 
lossy  transmission  lines  to  represent  material  damping  in  the  AEF  scheme  and  silver  epoxy 
bondlines.  The  effect  of  these  additions  was  to  significantly  dampen  ultrasonic  standing 
waves  created  in  the  AEF  system  which  resulted  in  reduced  power  transfer  efficiencies 
compared  with  theoretical  results  reported  intheliterature.TheLTSpice model  was  used  to 
investigate  an  AEF  system  installed  on  aluminium  plate  and  the  simulation  results  were 
compared  with  the  results  from  the  experi  ment.  M  odel  I  i  ng  and  experi  mentati  on  was  carri  ed 
outwithtwo  piezoelectric  elements  (transmit  and  receive)  thatwerewell  matched,  meaning 
thatthey  had  near  identical  resonance  and  anti -resonant  thickness  modes.  Itwasshownthatif 
the  drive  frequency  was  equal  to  the  anti-resonant  frequency  of  the  piezoelectric  elements 
used  to  create  the  AEF  system,  then  the  power  transfer  efficiency  was  unaffected  by  plate 
thickness.  Using  Pz27  disks,  of  38  mm  diameter  and  2  mm  thick,  laboratory  testing  has  shown 
thatwith  1W  of  input  power  then  300mW  of  real  power  could  be  transmitted  to  an  electrical 
load  through  aluminium  plates  of  1.6mm  and  5mm  thickness.  Futureworkwill  involveboth 
the  i  nvesti  gati  on  of  vari  ous  types  of  pi  ezoel  ectri  c  materi  al  s  to  determi  netheopti  mum  cerami  c 
composition  for  power  transfer  and  the  development  of  AEF-based  communications 
techniques. 


23 


DSTO-RR-0338 


5.  References 


1.  S.  C.  Galea,  S.  van  der  Velden,  S.  Moss,  I.  Powlesland,  "On  the  way  to  autonomy: 
the  wireless-interrogated  and  self-powered  'smart  patch' system,"  (to  be  published  in 
Encyclopaedia  of  Structural  Health  Monitoring  by  Wiley) 

2.  S.  D.  Moss,  S.  C.  Galea,  I.  G.  Powlesland,  M.  Konak,A.  A.  Baker,  "In-situ  health 
monitori  ng  of  a  bonded  composite  patch  usi  ng  thestrai  n  ratio  technique,”  SPI  E's  2000 
Symposi  u m  on  Smart  M  ateri al s  and  M  EM  S,  Smart  Structu res  and  Devi ces  Conference, 
Melbourne,  Australia,  SPIE  4325,  2000;  Paper  4235-41. 

3.  S.  C.  Galea,  S.  D.  Moss,  I.G.  Powlesland,  A.  A.  Baker,  "A  ppli  cation  of  a  Smart  Patch 
on  an  F/A-18  Aileron  Hinge,"  Proc  ACUN  -  4  Composite  Systems  - 
Macrocomposites,  Microcomposites,  Nanocomposites,  UNSW,  SYDNEY,  Australia, 
2002. 

4.  C.  Davis,  W.  Baker,  S.  M oss,  S.  Galea,  R. Jones,  "In  situ  health  monitoring  of  bonded 
composite  repairs  using  a  novel  fiber  Bragg  grating  sensing  arrangement," 
Proceedings  of  SPIE  --  Volume  4934,  Smart  Materials  II,  Alan  R.  Wilson,  Editor, 
November  2002,  pp.  140-149 

5.  N.  Rajic,  "In  situ  Monitoring  of  Crack  Growth  in  Mild  Steel  under  Closure 
Conditions  Using  a  Piezotransducer  Array,"  Journal  of  Intelligent  Material  Systems 
and  Structures,  Vol.  11,  No.  9,  696-702  (2000) 

6.  R.  Bartholomeusz,  B.  Crawford,  C.  Davis,  S.  Galea,  B.  Hinton,  G.  McAdam,  I. 
Powlesland,  K.  Sharp,  P.  Trathen,  T.  Trueman  and  A.  Wilson,  "DSTO  Aircraft 
Structural  Prognostic  Health  Monitoring  Program  for  Corrosion  Prevention  and 
Control,"  Corrosion  &  Prevention  -  2004  Australasian  Corrosion  Association  Inc.,  Perth, 
N  ovember  21-24  2004,  Perth,  WA 

7.  C.  P.  Gardiner,  P.  Vincent,  A.  Wilson,  D.  Ellery,  and  T.  Armstrong,  "A  trial  sensor 
network  for  the  Armidale  Class  Patrol  Boat,”  Pacific  2008  International  Maritime 
Conference,  Sydney,  Australia,  January  2008 

8.  P.  S.  Vincent,  P.  J.  McMahon,  R.  F.  Muscat,  L.  Zeve,  A.  R.  Wilson,  "A  small  low- 
power  networked  and  versatile  sensor  interface  (Invited  Paper),"  SPIE  International 
Symposium  on  Smart  Materials,  Nano-,  and  Micro-Smart  Systems  2006:  Smart  Structures, 
Devices  and  Systems  III,  December  2006,  [6414-39] 


24 


DSTO-RR-0338 


9.  Y.  Hu,  X.  Zhang,  J.  Yang,  Q.  Jiang,  "Transmitting  Electric  Energy  Through  a  Metal 
Wall  by  Acoustic  Waves  Using  Piezoelectric  Transducers,"  IEEE  Trans.  Ultrason., 
Ferroelect.,  Freq.  Cont.,  vol.  50,  no.  7,  pp.  773-781,  2003 

10.  S.  Sherrit,  M .  Badescu,  X.  Bao,  Y.  Bar-Cohen,  Z.  Chang,  "Efficient  electromechanical 
Network  Models  for  Wireless  Acoustic- Electric  Feedthroughs,"  Proceedings  of  SPIE 
Smart  Structures  Conference,  San  Diego,  CA.,  Mar  6-10,  2005,  SPIE  vol.  5758 

11.  M.  Engelhardt,  LTspice'  SwitcherCAD  III,  Linear  Technology  Corporation, 
www.linear.com ,  2008 

12.  S.  A.  Morris,  C.  G.  Hutchens,  "Implementation  of  Mason's  Model  on  Circuit 
Analysis  Programs,"  IEEE  Trans.  Ultrason.,  Ferroelect.,  Freq.  Cont.,  vol.  U  FFC-33,  no.  3, 
pp.  295-298, 1986 

13.  W.  M.  Leach  Jr.,  "Control  led-Source  Analogous  circuits  and  SPICE  Modelsfor 
PiezoelectricTransducers,"  IEEE  Trans.  Ultrason.,  Ferroelect.,  Freq.  Cont.,  vol.  41,  no.  1, 
pp.  60-66, 1994 

14.  A.  Puttmer,  P.  Hauptmann,  R.  Lucklum,  0.  Krause,  B.  Henning,  "SPICE  Model  for 
Lossy  PiezoceramicT ransducers,"  IEEE  Trans.  Ultrason.,  Ferroelect.,  Freq.  Cont.,  vol .  44, 
no.  1,  pp.  60-66, 1997 

15.  J .  van  Deventer, T  Lofqvist,  "PSpice Simulation  of  U Itrasonic Systems,"  IEEE  Trans. 
Ultrason.,  Ferroelect.,  Freq.  Cont.,  vol.  47,  no.  4,  pp.  1014-1024,  2000 

16.  L.  Wu,  Y-C.  Chen,  "PSPICE  approach  for  designing  the  ultrasonic  piezoelectric 
transducer  for  medical  diagnostic  applications,"  Sensors  and  Actuators  A,  vol.  75,  pp. 
186-198, 1999 

17.  W.  P.  M  ason  Electromechanical  Transducers  and  Wave  Filters,  lstedn  (Princeton,  NJ: 
Van  Nostrand)  1948 

18.  R.  Krimholtz,  D.  A.  Leedom  and  G.  L.  Matthaei,  "New  equivalent  circuit  for 
elementary  piezoelectric  transducers,"  Electron  Lett,  .vol.  6,  pp.  398-9, 1970 

19.  M.  Kim,J.  Kim,  W.  Cao,  "Aspect  ratio  dependence  of  electromechanical  coupling 
coeff i  ci  ent  of  pi  ezoel  ectri  c  resonators, "  App.  Phys.  Lett. ,  vol .  87,  pp  132901- 1  to  132901- 
3,  2005 


25 


DSTO-RR-0338 


20.  D.  A.  Berlincourt,  D.  R.  Curran  and  H.  Jaffee,  Physical  Acoustics  Principles  and 
Methods,  Volume  I- Part  A,  pp. 170-267  (Academic  Press,  New  York  and  London,  1964), 
Library  of  Congress  Catalog  Card  N  umber:  63-22327 

21.  Private  communication,!.  M  adsen,  Ferroperm  Quality  M  anager,  4th  October,  2007 

22.  Piezoelectric  material  parameters,  Ferroperm  website,  www.ferroperm.com, 
September  2006 

23.  C i  rcu i tw or ks  Conductive  Epoxy,  CW2400,  http:/  /  www.chemtronics.com,  2008 

24.  N .  Rajic,  "A  model  forthepiezoelectrictransduction  of  stress  waves,”  Smart  Mater. 
Struct.,  vol.  15,  pp.  1151-1164, 2006 

25.  http:/  /  www.fairchildsemi.com/  pf/  IN/  lN5818.html,  2008 

26.  W.  P.  M  ason,  Physical  Acoustics  Principles  and  Methods,  Volume  I -Part  A,  pp. 485-499 
(Academic  Press,  New  York  and  London,  1964),  Library  of  Congress  Catalog  Card 
N  umber:  63-22327 

27.  J.  Johansson,  P.  Marti  nsson,  "Incorporation  of  Diffraction  Effects  in  Simulations  of 
Ultrasonic  Systems  using  Pspice  models,"  Proc.  2001  IEEE  Ultrasonics  Symposium,  vol. 
1,  pp.  405-410,  Atlanta,  GA,  USA,  7-10 Oct.  2001 

28.  G.  S.  Kino,  Acoustic  Waves:  Devices,  Imaging  and  Analog  Signal  Processing, 

Englewoods  Cliffs,  NJ:  Prentice- Hal  1, 1988 

29.  H.  Seki,  A.  Granato,  R.  Truell,  "Diffraction  Effects  in  the  Ultrasonic  Field  of  a 
Piston  Source  and  Their  Importance  in  the  Accurate  Measurement  of  Attenuation," 
The  J.  Acoust.  Soc.  Amer.,  vol.  28,  no.  2,  pp230-238,  M  ar.  1956 


26 


DSTO-RR-0338 


Appendix  A:  Example  piezo  parameter  calculation 


To  calculate  f33s:  given  that  the  permittivity  e33T  =s33,rT  e0  van  Deventer  et  al  [15]  gives  the 
following  relationship  between  e33J  and  e33s, 


4  1 

4  1-4J2 


(Al) 


Because k33  and  e33T  are  normally  provided  by  the  manufacturer,  e33s  may  calculated  using 
equation  A 1. 

To  calculate  e33:  piezoelectric  constant  e33  may  be  calculated  using, 

c33  =  2  d^ i  c13  +  d 33  c33  (A 2) 


Equation  A2can  beused  ifthe  manufacturer  specifies  c23£and  c33E.  If  these  elastic  stiffness's  are 
not  provided  then  e33  may  be  estimated  using, 


e33  —  r/33  F33  (kT ) 


(A3) 


Equations  A2  and  A3  can  be  used  to  calculate  an  accuratee33  for  most  types  of  PZT  ceramic, 
however  the  calculated  e33  can  be  erroneous  for  PZT  with  unusual  microstructure  (e.g.  Pz21 
and  Pz46). 


To  calculate  c33D:  The  relationship  between  c33°  and  c33E  is, 


(A  4) 


A  PZT  manufacturer  normally  specifies  both  kT  and  c33E  allowing  c33D  to  be  found.  If  the 
manufacturer  does  not  specify  c33E  then  it  will  need  to  be  found  from  Y33E, 


c 


E 

33 


1  -A2 


(A  5) 


Assumi  ng  that  kT,  k33,  and  kP  are  known  for  a  material,  the  constant  A  can  be  found  from  the 
relationship  given  by  Berlincourtef  al  [20], 


kT 


ki3  -  Akp 

Vi-a27i-4 


A  can  also  be  found  from  its  definition, 


(A  6) 
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A  = 


V2  sf3 

V4(4  +4) 


(A  7) 


wheres23E,  s33E,  suE  and  suE  are  compliances  at  constant  field. 

With  C33°  known  then  vD,  the  clamped  speed  of  sound  in  the  piezoelectric  material,  can  be 
found  using, 


(A  8) 


Also  required  for  lumped  element  modelling  are  the  density  and  speed  of  sound  of  the 
medium  that  the  PZT  element  is  situated  in.  For  the  example  shown  in  Figure  la  the 
surrounding  medium  is  air,  which  has  density  p  =  1.25  kg/  m3  and  speed  of  sound  v 
=  343.4  ml  s .  K  nowing  the  area  AE  of  the  PZT  element  allows  the  mechanical  loading  dueto 
the  medium  to  be  found, 

Z0=pvAE  (A9) 

Using  a  single  resistor  Z0  to  model  the  medium  is  equivalent  to  making  the  assumption  that 
no  acoustic  echo  will  be  received  back  from  the  medium  [16], 


28 


DSTO-RR-0338 


Appendix  B:  LTSpice  netlists 


Given  below  are  example  of  LTSpice  library  and  schematic  files  used  to  perform  the 
modelling  shown  inthispaper.  N  ote  that  the  parameters  given  in  PZTtrans.libaredefined  in 
Appendix  C. 

B.l.  PZTtrans.lib 


*PUTTMER  1997 THICKNESSMODE  TRANSDUCER 

.SUBCKT  PZTtrans  E  B  F 

01B1F  lmyLossyTlineTRANS 

VI 1  2 

El  2  040  1 

V2E  3 

CO  3  0  {cOT} 

FI  0  3  VI  fhcOT} 

F204V2  -{hT} 

R1401E3 

C1401 

.model  myLossyTlineTRANS  LTRA(len=-8;hickT}R={rT}L=-!LT}C={cT}) 
.ENDS 

B.2.  PZTtrans.asy 


Version  4 

SymbolType  BLOCK 
RECTANGLE  Normal  48  48-48-48 
TEXT  -48  -1  Left  0  PZTtrans 
SYMATTR  Value  PZTtrans 
SYMATTR  Prefix  X 
SYMATTR  SpiceModel  PZTtrans.lib 
SYMATTR  Value2  PZTtrans 
PIN  -48-32  LEFT  8 
PINATTR  PinNameE 
PINATTR  SpiceOrder  1 
PIN  -48  32  LEFT  8 
PINATTR  PinNameB 
PINATTR  SpiceOrder  2 
PIN  48  32  RIGHT  8 
PINATTR  PinNameF 
PINATTR  SpiceOrder  3 
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B.3.  const_power.lib 

.subckt  const_power  in  m_Power 
B1  N0010V=V(vd)*sin(2*pi*ti  m^V  (f )) 

R1  in  N 001 50 

*for  real  power  the  next  I  i  ne  must  be: 

*B2  zlOV— l*V(in)*i(rl) 

*for  apparent  power  the  next  I  i  ne  must  be: 
*B2zlOV=abs(V(in)*i(rl)) 

B2  zl  0  V  =abs(V  (i  n)*i  (rl)) 

R3  m_Power  zl  IK 
C2  m_Power  0  20n 
Vlf  OfDRIVEcp} 

V2  Gain  0100 

V3  ldeal_Power  0  -{pDRIVEcp} 

B4  0vd  l=V(Gain)*(V(ldeal_Power)-V(m_power)) 
R4  vd  0  le8 
C3vd  OlOp 
.icV(vd)=10 

.ic  V(mjDower)=-{pGUESScp} 

.backanno 

.ends 
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Appendix  C:  Parameter  values  used  in  calculations 


definition 

designation 

units 

material 

Hu  PZT 
[9] 

Pz27 

[22] 

aluminium 

[5] 

steel 

[9] 

silver 
epoxy  [5] 

Clamped  relative 
permittivity 

£S 

°33  ,rel 

- 

1470.5 

914 

- 

- 

- 

piezoelectric 

constant 

e33 

C/  m2 

23.3 

16 

- 

- 

- 

clamped  elastic 
modulus 

c33 

GPa 

159 

144 

73.1 

269 

7.28 

density 

P 

kg/  m3 

7500 

7700 

2770 

7850 

3890 

mechanical 

quality 

Q 

- 

136 

81.5 

1000 

100 

23.4 

Definition 

Designation 

Equation 

Effective  Units 

permittivity  of  free  space 

e0 

8.8542*1(1 12 

C/  (V.m) 

clamped  permittivity 

eST 

e33S*eO 

C/  (V.m) 

piezoelectric  constant 

e33T 

C/  m2 

clamped  elastic  modulus 

C33DT 

N/  m2 

density 

rhoT 

kg/  m2 

mechanical  Q 

QT 

- 

disk  thickness,  and  radius 

thickT,  radT 

m 

clamped  velocity  of  sound 

vDT 

SQRT(c33DT/  rhoT) 

m /  s 

anti -resonant  frequency 

resfreqT 

vDT/  (2*thickT) 

Hz 

disk  area 

areaT 

pi*radT*radT 

m2 

capacitance 

cCTT 

eST*areaT/  thickT 

F 

transmitting  constant 

hT 

hear 

e33T/  eST 
hT*COT 

V/  m 

F.V/  m 

acoustic  imped ance*area 

zcrr 

rhoT*VDT*areaT 

kg/  s 

lossy  transmission  line  inductance 

LT 

z  CT/  vDT 

kg/  m 

lossy  transmission  line  capacitance 

cT 

V  (vDT*zOT) 

s/  kg 

lossy  transmission  line  resistance 

rT 

2*p  i  *resf  reqT*LT /  QT 

kg/  (m.s) 

density  of  medium  back  face 

rhoTB 

kg/  m3 

velocity  of  sound  in  back  medium 

vDTB 

m /  s 

acoustic  impedance*area 

zOTB 

rhoTB*vTB*areaT 

kg/  s 

density  of  medium  front  face 

rhoRB 

kg/  m3 

vel  oci  ty  of  sou  nd  i  n  front  med  i  u  m 

vDRB 

m/  s 

acoustic  impedance*area 

zORB 

rhoRB*vRB*areaT 

kg/  s 
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